On the origin of the particles in black hole evaporation 
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We present an analytic derivation of Hawking radiation for an arbitrary (spatial) dispersion re- 
lation uj{k) as a model for ultra-high energy deviations from general covariance. It turns out that 
the Hawking temperature is proportional to the product of the group dcu/dk and phase u/k veloci- 
ties evaluated at the frequency u of the outgoing radiation far away, which suggests that Hawking 
radiation is basically a low-energy phenomenon. Nevertheless, a group velocity growing too fast at 
ultra-short distances would generate Hawking radiation at ultra-high energies ("ultra-violet catas- 
trophe" ) and hence should not be a realistic model for the microscopic structure of quantum gravity. 
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Introduction In the history of physics, unexplained 
coincidences were often the precursor of striking discov- 
eries. For example, the limiting propagation velocity de- 
rived from the properties of coils and capacitors in table- 
top experiments turned out to be strikingly close to the 
speed of light measured via planetary motion - which 
guided the unification of these phenomena by the laws 
of electrodynamics and ultimately lead to the theory of 
relativity. Today, a similar riddle is the question of why 
black holes seem to behave like thermal objects [1| and 
evaporate by emitting Hawking radiation y, y, IJ] . 

One way of achieving a better understanding of these 
links is to study the origin of the particles in black hole 
evaporation, i.e., the question of where they are created. 
For example, it has been suggested to resolve the black 
hole "information paradox" (i.e., the apparent contradic- 
tion between unitarity and the second law of thermody- 
namics in this system) by encoding information into the 
outgoing Hawking particles. Clearly, this (hidden) en- 
coding mechanism should then occur at the origin of the 
particles (or on their way to infinity). 

Since the event horizon marks the "point of no re- 
turn" , the modes containing the Hawking particles orig- 
inate from the region very close to the horizon, i.e., from 
very short wavelengths (gravitational red-shift). How- 
ever, the origin of the modes is not necessarily the place 
where the particles are created. 

In order to address this problem, we derive Hawking 
radiation in the presence of a very general dispersion re- 
lation uj = ck —f CO = uj(k) associated to the propa- 
gating degrees of freedom. Such a modified dispersion 
relation is inspired by the analogy to condensed mat- 
ter, i.e., the sonic black hole analogues (silent or "dumb" 
holes), which rely on the quantitative analogy between 
quantum fields in curved space-times and phonons (or 
other quasi-particles) propagating in fluids with a gen- 
eral flow velocity v{t, r), cf. [5, 6]. In this case, the phase 
uj/k and group duj/dk velocities vary with wavelength 
and thus the dependence of Hawking radiation on the 



dispersion relation should show us which wavenumbers 
k are most important for particle creation: The emit- 
ted radiation at various frequencies to will "see" different 
horizons. The question this paper will try to answer is 
what determines the temperature of the radiation emit- 
ted at any particular frequency uj. If the temperature 
is determined when the wavelengths are very small, and 
frequencies large - i.e., when the horizon first splits the 
incoming wave packet (in its vacuum state) into posi- 
tive pseudo-norm modes outside the horizon (which will 
turn into the Hawking particles) and negative pseudo- 
norm modes (their infalling partner particles) inside - 
then one would expect a universal temperature for all 
of the low-frequency modes. On the other hand, if it is 
the low frequency aspects of the modes which determine 
the temperature, one might expect the properties of the 
horizon deflned via one of the velocities (phase, group, or 
other) associated with the wave at low frequencies to dic- 
tate the temperature (via dv/dr at that horizon location) 
of that mode. 

Previous analytic calculations were restricted to low 
energies to and a small vicinity of horizon, see, e.g., 
[7|, while numerical studies necessarily involved a given 
(mostly sub-luminal) dispersion relation within a re- 
stricted parameter range. In the following, we present 
an analytic derivation of Hawking radiation valid for any 
frequency lo and almost arbitrary (spatial) dispersion re- 
lations uj{k); the only assumption we make is that the 
black or "dumb" hole should be large, i.e., macroscopic, 
and that the velocity proflle of the background flow have 
a specific form. 

Dispersion relation Since Hawking/ "dumb" -hole ra- 
diation is basically a 1+1 dimensional effect, we consider 
the Painleve-GuUstrand-Lemaitre metric in 1+1 dimen- 
sions {h~c~ Gn = fee = 1) 



ds"^ = [l - u^(.t)] dt^ - 2v{x) dtdx ~~ dx^ 



(1) 



with the local frame dragging velocity v(x)^ which cor- 
responds to the flow velocity of the fluid analogue [51]. 



It determines the position of the horizon via v{x) = ±c 
where c is the velocity of the propagating modes, which 
is assumed to be constant and set to unity here. The 
propagation of a massless scalar field <& in this metric 
(or, equivalently, phonons in the fluid) is governed by 
the d'Alembertian 
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(2) 



Due to the conformal invariance in 1+1 dimensions, the 
left-moving modes {dt + d^v — d^)^ = propagating 
against the frame-dragging velocity are decoupled from 
the right-moving solutions (p = {dt + d^v — dx)^- How- 
ever, an arbitrary modification of the dispersion relation 
would not preserve this decoupling in general. The choice 



ah^{dt + dxv){dt + vdx)-h{dl). 



(3) 



for example, does not factorize for a general function h. 
On the other hand, if we modify the left and right- moving 
branch separately with an arbitrary function / via 
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f = {dt + dx[i + v + fi~dl)])x 

{dt-[l-v + f{-dl)]dx) , 



(4) 



the left-moving modes {dt— [l — v + f{—dx)]dx)^ = are 
again decoupled from the right-moving solutions given by 
0= {dt-{\-v^j{~dl)\dx)^. The difference between the 
two options ([3]) and ^ scales with the derivative of v{x) 
compared with the characteristic scale of the dispersion 
relation \^{x),j{—d'^\, which is negligibly small in the 
hydrodynamic limit (which corresponds to the black hole 
being large). The coupling occurs at high frequencies 
over long scales (the scale of variation of v) which means 
the response will be adiabatic and left-moving waves will 
not convert to right-moving. This is the reason why the 
numerical simulations using the first option Q never saw 
a mixing between left and right-movers one might expect 
from the coupling between the two sets of modes. 

Derivation Even though the Schwarzschild metric 
corresponds to v{x) — ±^J2M /x, we consider two differ- 
ent velocity profiles v{x) = —X/x and v{x) = kx in the 
following, because they admit analytic solutions. The 
latter of course has no asymptotically flat region, and 
calling it a black, or "dumb", hole is metaphorical. Let 
us first study the case v{x) = —X/x. After a Fourier- 
Laplace transformation with dx — > ik and x —^ idk as 
well as dt -^ —iuj, the left-moving solutions satisfy the 
integral equation (since d^^ = f dk) 
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/(fc2)])0^(fc) = O, 



(5) 



which can be solved via separation of variables and dif- 
ferentiation (l)u){k) = exp{—iX J dk' / g{k')} / g{k) with the 
spectral function g{k) = 1 + /(fc^) — uj/k. The inverse 
Fourier-Laplace transformation 



^Ur) = 



dk 



exp < ikx — iX 



dk' 



(6) 



yields the spatial modes (j)uj{x) with the integration con- 
tour being determined by the boundary conditions. In 
the following, we shall assume the length scale A on which 
v{x) changes to be very large compared with the typical 
wavenumber of the dispersion relation f{k^) and the fre- 
quency cu. In terms of the fluid analogue, this is precisely 
the hydrodynamic limit - whereas, for real black holes, is 
corresponds to demanding that the size of the black hole 
is much larger than the Planck scale. Note, however, that 
we do not restrict ui relative to the Planck scale (i.e., cu 
could be Planckian). Since the exponent in the integral 
above contains the large numbers x and A, it is very use- 
ful to deform the integration contour into the complex 
plane (assuming that / is an analytic function), where 
the leading contributions will be determined by singu- 
larities and the associated branch cuts as well as saddle 
points. The saddle points (stationary phase) 
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f{k' 



(7) 



are solutions of dispersion relation oj + vk = k[l + /(A;^)]. 
For large wavenumbers |fc| 3> w, we get pairs of saddle 
points k± satisfying /(fc^) = X/x — 1. The singularities 
of the integrand at g{k) = correspond to solutions of 
dispersion relation far away from the black hole a; — » cxd. 
Assuming simple poles at ka only, we may employ the 
residual expansion 
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I f dk ka 

2ni J g{k) Vg,.{ka) 



(8) 



The contours £q denote small circles around the poles 
at ka and the residual coefficients Cq are related to the 
group velocity at these points. Insertion of the residual 
expansion ([8]) into Eq. ^ yields 



<^"W = /Jy-^'^n(^-M' 



(9) 



Consequently, there are branch cuts starting from each 
singularity unless 2iMca G Z, cf. Fig. [1] 
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FIG. 1: [Color online] Sketch (not to scale) of integration 
contours in the complex plane for the sub-luminal (left) and 
the super- luminal case (right) . The dots denote saddle points 
k± at the real axis and the squares are the singularities ka 
with the associated branch cuts (dashed lines). 



Sub-luminal case In order to determine the most suit- 
able integration contour in the complex plane, we have to 
incorporate some assumptions about the function f{k^)- 
First, we assume that the asymptotic (a; — > cxd) disper- 
sion relation uj = fc[l -I- /(fc^)] is convex, i.e., sub-luminal 
/(fc^) < 0, and always monotonically increasing, i.e., 
with a positive group velocity. In this case, saddle points 
in Eq. ([7|) occur outside the horizon x > X and yield the 
following contributions to the integral in Eq. ([6]) 



'/'^(a^) 
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Xg'{k±) 



e^'=±^ [] (fc± - fc„) 



-iXCa 



(10) 



In view of their spatial behavior, these are the positive 
(j)'^ and negative 0J pseudo-norm in-modes with large 
wavenumbers \k±\ 3> w and therefore small group veloc- 
ities, which are swept towards the horizon. After quan- 
tizing the field <&, these positive/negative pseudo-norm 
solutions yield the creation/annihilation operators of the 
in-modes at large k. 

In order to close the integration contour, we have to 
circumvent the branch cuts in the upper complex half 
plane 3(fc) > 0, cf. Fig. [1] In the limit x -^ oo, the 
contributions of the branch cuts starting from the singu- 
larities ka away from the real axis '^{ka) > are expo- 
nentially suppressed and only the branch cut starting at 
the real axis 9(fc*) = contributes. This wavenumber 
fc* represents a real solution of the dispersion relation far 
away from the horizon and just corresponds to the out- 
going Hawking radiation with frequency cj > 0. Hence, 
at cc — > oo, the integral in Eq. ^ corresponding to the 
contour €y in Fig. [T] yields a superposition of the outgo- 
ing Hawking modes with fc* and the large-fc in-modes in 
(fTO| with k±. Continuing this solution beyond the hori- 
zon X < A, the saddle points vanish and the integrand in 
()10p decays exponentially in the lower complex half plane 
3(fc) < 0. Thus, we deform the integration contour to 
£<, where the main contribution stems from the branch 
cut(s) starting at ^{ka) < 0. Again, for large A, these 
contributions are exponentially suppressed. 

Therefore, the outgoing Hawking mode with fc, origi- 
nates entirely from the in-modes in (jlOp . Assuming that 
the initial quantum state is the ground state of the large- 
k modes in (fTO)) with respect to the freely falling observer, 
the amount of created particles is then determined by 
the mixture between these positive and negative pseudo- 
norm solutions with large k contained in the outgoing 
fc*-modes after the immense gravitational red-shift near 
the horizon. In view of \ka\ ^ \k±\, the only differ- 
ence between positive and negative pseudo-norm modes 
is caused by branch cut(s) 
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exp<^7rA^(-l)""5R(ca) 



(11) 



where (—1)*° is the sign associated with the direction of 
the branch cut. Since g{k) is a real function, the singu- 



larities ka occur symmetric w.r.t. the real axis ka -^ k^ 
and Cq -^ c*. Choosing the branch cuts suitably (see 
Fig. [T]) , the contributions from the symmetric pairs can- 
cel each other and hence only the singularity at the real 
axis fc* = fed G K contributes. Together with the uni- 
tarity relation la^^p — |/3cjp = 1, the ratio pTjl directly 
determines the size of the Bogoliubov coefficients via 
\auj/Pio\'^ = \(l>i/4>Z\'^ = exp{w/r}. Hence we may read 
off the effective Hawking temperature 



2Hawking(^) — 



fgr(fc*)^ph(fc*) 

2ttX 



(12) 



We observe that the geometric mean of group and phase 
velocity [8| evaluated at the frequency uj of the outgo- 
ing radiation far away x — > oo determines the Hawking 
temperature [9|. Therefore, the behavior of the disper- 
sion relation at large k is not relevant - even though the 
Hawking radiation originates from large-fc modes - which 
indicates that the Hawking effect is basically a low-energy 
phenomenon. The cj-dependence of the Hawking tem- 
perature can be explained by the fact that high-energy 
wave-packets have a different group velocity than those 
at low energy and hence the various modes "see" differ- 
ent horizons and thus other values for the surface gravity, 
i.e., velocity gradient dv/dx = Xjx^ (x v^/X. 

This explanation has been confirmed by numerical sim- 
ulations [9] and can further be supported by considering 
the second case v{x) = nx, where the velocity gradient n 
is constant. In this case, Eq. ^ should be replaced by 

(t>u:{x)= f dkk-"^/''exp\ikx-iK fdk'f{k'^)\ .(13) 



Hence the weight of the branch cut starting at fc = is 
just determined by the ratio lu/k and the Hawking tem- 
perature does not depend on to at all because all modes 
"see" the same surface gravity k 



Tkawking — 7;— — COUSt . 

Ztt 



(14) 



Super-luminal case A super-luminal dispersion rela- 
tion f{k'^) can be treated in a completely analogous way. 
As the only difference, the large-fc in-modes determined 
by the saddle points k± originate from inside the horizon 
X < X and thus the contours in the complex plane (£> for 
X > X and £< for a; < A, cf. Fig. [1]) are slightly different. 

In view of this observation, one might wonder whether 
the ansatz v{x) = —X/x instead of v{x) = ± y^2M/x 
is justified. To address this question, let us consider 
the Schwarzschild geometry in 1-f 1 dimensions using the 
Eddington-Finkelstein coordinates (V, r) 



ds' 
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2dVdr. 



(15) 



For a massless scalar field $, the wave equation reads 



2dvdr + dr 



r 



a,. $ = , (16) 



where we have again included a modification /(fc^) of 
the dispersion relation. Comparison with the previous 
derivation yields completely the same results for the out- 
going solutions (p = dr^ wp to the replacement u ^ 2uj 
due to the Eddington-Finkelstein coordinate t ^ V. 

However, in the super-luminal case, an additional com- 
plication may arise: According to Eqs. ([8]) and pTj) . 
the thermal Boltzmann factor exp{a;/(4T)} determining 
the amount of created particles with frequency lu can 
be recast into the alternative form exp{27rMfc*/i'gr(A:*)}. 
Hence, if the group velocity grows slower than linear in 
k, the number of produced particles decreases with en- 
ergy. However, if Wgi(^) rises too fast in some fc- region, 
the amount of created particles drops at low k (where 
fgi- = 1) but later increases again! In such an extremal 
case, the Hawking radiation could contain a large contri- 
bution of ultra- high energy particles ( "ultra-violet catas- 
trophe" ) . Going a step further and taking the dispersion 
relation seriously as a model for ultra-high energy de- 
viations from general relativity [lO|, one would exclude 
such a case in view of our observational evidence for the 
existence of black holes with macroscopic life-times. 

Let us discuss some examples: The dispersion rela- 
tion oj'^ — k'^ + fc^ which is realized for the sonic black- 
hole analogues in Bose-Einstein condensates 11|, does 
not generate an "ultra-violet catastrophe" and repro- 
duces Hawking's prediction. In contrast, the expressions 
w^ = exp{fc^} — 1 or cj = k/^/Y^^W grow too fast and 
hence lead to the aforementioned problems [12| . 

Note that the condition Vgr{k*) > 0{Mk^) for par- 
ticle creation obtained from the Boltzmann factor pre- 
cisely marks the break-down of the saddle-point (i.e., ge- 
ometric optics) approximation. Writing the integrand in 
Eq. ([6]) as G'(fc) exp{i^(fc)}, the first-order corrections to 
the saddle-point expansion scale as G' /G x _F("^)/(_F")^, 
i^(4)/(^/')2^ and G"/{GF"), evaluated at the saddle 
point F' = 0. In our case (O, we have F" = -2iMG' 
and hence inserting Vg,:{k^,) > 0{Mk^,) yields "correc- 
tions" of order one - i.e., the saddle-point approximation 
fails. Even if the modes started out in their ground state 
(at k±), they get excited (at k^) due to the strong gravi- 
tational red-shift. Based on these general adiabaticity ar- 
guments, one would expect that the main result remains 
valid even beyond a sole modification of the dispersion 
relation: If the outgoing Hawking modes originate from 
the vicinity of the singularity and the spectral proper- 
ties of quantum gravity change too fast with energy, one 
would expect a break-down of adiabaticity at short dis- 
tances resulting in the emission of high-energy particles. 
This mechanism is not necessarily restricted to a Planck- 
length vicinity of the singularity since the effective sur- 
face gravity "seen" by the high-energy modes scales with 
2M/r^ and hence may exceed the Planck temperature 
already many Planck lengths away from the singularity. 
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In order to motivate the third example, one could ar- 
gue that it should be impossible to resolve spatial de- 
tails smaller than the Planck length £p = ^JhGti/c? ~ 
1.6 X 10~^^m, because Heisenberg's uncertainty prin- 
ciple would then require concentrating a large (trans- 
Plankian) energy and thus mass in a small volume - 
which should form a black hole. In analogy to special rela- 
tivity, where the acceleration of massive particles towards 
the speed of light is prohibited by a diverging energy 
E = m<? I ^\ — v^ l(? , one could speculate that quan- 
tum gravity impedes a spatial resolution better than the 
Planck length i-p in a similar manner lo — ck/ yjl — £pk^. 
Even though the two relations cj^ — k'^ + k and uj = 
k/\/l — fc^ can probably not be distinguished via obser- 
vations in the low-energy sector (e.g., energy-dependent 
time-delay measurements of signals from distant galax- 
ies), they induce drastic differences for black holes. 



